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Abstract 



In this paper we provide an explicit construction of any arbitrary unitary transformation on n qubits 
from one qubit and a single two qubit quantum gate. Building on the previous work demonstrating 
the universality of two qubit quantum gates, we present here an explicit implementation. The con- 
struction is based on the Cartan decomposition of the semi-simple Lie group SU(2 n ) and uses the 
geometric structure of the Riemannian Symmetric Space su(2™- i )®su(2 n - i )®u(i) ■ decomposi- 
tion highlights the geometric aspects of the problem of building an arbitrary unitary transformations 
out of quantum gates and makes explicit the choice of pulse sequences for the implementation of 
arbitrary unitary transformation on n coupled | spins in NMR quantum computing. Finally we 
make observations on the optimality of the design procedure. 

a 
cr. 



1 Introduction 



Recent interests in quantum computing have been fostered by the prospects of building a compu- 
tational theory which is different and more powerful than the classical computational theory. The 
vision of solving a class of problems using quantum computing, which have been known to be com- 
putationally hard in the classical setting, has been advanced by the recent work of |jl3], Any 
quantum computation involves evolution of the initial state under a series of unitary transformations. 
Absolutely essential to the design of an universal quantum computer is the ability to synthesize any 
arbitrary unitary transformations from simple quantum gates (lower dimensional unitary transfor- 
mations). The search for such universal components was initiated by Deutsch's original discovery 
p0| of a three-bit universal quantum logic gate. Following which DiVincenzo pjj showed that two- 
bit universal quantum gates also exist. This result was then generalized by the work of Barenco et al. 
PH , Lloyd Q and Deutsch et al. |^0) showing that almost any two bit gate is also universal. The 
central idea behind demonstrating the universality of these gates rests on the fact that the genera- 
tors of the unitary transformations produced by these gates generate the Lie algebra of the unitary 
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group and therefore by suitable composition of these gates, we can produce any arbitrary unitary 
transformation. This is basically the notion of controllability in mathematical control theory- that 
is, whether available Hamiltonians can prepare an arbitrary state of the quantum system |24|, p5| . 
The conditions of controllability of a quantum system were rediscovered as the conditions for the 
universality of a quantum computer. 

Problems of similar nature arise in coherent spectroscopy. Many areas of spectroscopic fields, such 
as nuclear magnetic resonance (NMR), electron magnetic resonance and optical spectroscopy rely 
on a limited set of control variables in order to create desired unitary transformations |5|, |^, [?]]. 
In NMR, unitary transformations are used to manipulate an ensemble of nuclear spins, e.g. to 
transfer coherence between coupled spins in multidimensional NMR-experiments || or to implement 
quantum- logic gates in NMR quantum computers |@]. However, the design of a sequence of radio- 
frequency pulses that generate a desired unitary operator is not trivial ||. So far, no general 
design approach is known for the implementation of a desired unitary transformation During 
the last decade the questions of controllability of quantum systems have generated considerable 
interest |l7]]. In particular, coherence or polarization transfer in pulsed coherent spectroscopy 
has received lot of attention || Algorithms for determining bounds quantifying the maximum 
possible efficiency of transfer between non-Hermitian operators have been determined S. There is 
utmost need for design strategies for pulse sequences that can achieve these bounds. From a control 
theory perspective, all these are constructive controllability problems [ p~4[ . 

In this paper, we present a constructive solution to the problem of producing an arbitrary 2 n dimen- 
sional unitary transformation acting on n qubits from unitary transformations acting on individual 
and pair of qubits (elements of SU{2) and S*[/(4) respectively). The design strategy presented here 
is a generalization of well known Euler angle decomposition for rotations. Recall if 




represent the standard Pauli spin matrices, then any V 6 SU (2) has a decomposition 

V = exp(-ial x ) exp(-i/3/ 2 ) e-xp(-ijl x ). 
where a, j3, 7 <E R. Similarly any W € SU(2 n ), has a decomposition 

W = K 1 AK 2 . 

Where K 1 ,K 2 S SU{2 n ~ 1 ) <g> SU{2 n ~ 1 ) <g> U{1) and A G exp(J)), where f) is a Cartan subalgebra of 
the Riemannian symmetric space 

SU(2 n ) 

SU(2 n ~ 1 ) i SU(2 n - 1 ) ® U(l) ' 

We will elaborate on all these notions. The point to be emphasized here is that we then obtain a 
recursive formula for the decomposition. Given the decomposition W = K1AK2, we can further 
decompose K\ and K 2 in terms of the elements of SU{2 n ~ 2 ) <g) SU(2 n ~ 2 ) (g) U(l) and so on. Finally 
we will be left with only one and two qubit operations. We will see that this recursive formulation 
highlights the geometric structure of the problem. 
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We are also interested in the implementation of these unitary transformations in a network of 
coupled spins in the context of NMR. The emphasis being NMR quantum computers and coherence 
transfer experiments in multidimensional spectroscopy. We emphasize again that from a control 
theory viewpoint, these are constructive controllability problems. The dynamical system to be 
controlled is defined through the time-dependent Schrodinger equation 

U{t) = -iH(t)U(t), U(0)=I, 

where H(t) and U(t) are the Hamiltonian and the unitary displacement operators, respectively. In 
this paper, we will only be concerned with finite-dimensional quantum systems. In this case, we can 
choose a basis and think of H (t) as a Hermitian matrix. We can split the Hamiltonian 

m 
i=l 

where Hd is the part of Hamiltonian that is internal to the system and we call it the drift or free 
evolution Hamiltonian and Xa=i v i(t)Hi 1S the part of Hamiltonian that can be externally changed. 
It is called the control or rf Hamiltonian. The equation for U (t) dictates the evolution of the density 
matrix according to 

p(t) = U(t)p(0)U\t). 

where 

m 

U = -i(H d + J2viHi)U;U(0) = I . (1) 
i=i 

In a network of coupled spin i nuclei, the control terms Hi correspond to the Hamiltonian that effect 
each spin individually (we assume that resonance frequencies of the spins are well separated so that 
selective excitation is possible). The group K generated by the Lie algebra {Hi} la corresponds to 
n direct copies of SU(2), which we denote by 

[SU{2)f n . 

The drift term Hd is the part of the Hamiltonian that corresponds to couplings among the spins. 
It is known that if the Lie algebra {Hd, Hi\^A equals su(2"), the algebra of 2™ x 2™, traceless skew 
Hermitian matrices, then we can steer the system (Q), from identity to any point Uf G SU(2 n ). The 
problem we will address in this paper is explicit construction of a pulse sequence or in other words 
control laws u,*, that steer the system (|l|) from U(0) = I to some Uf, in finite time. 

We begin by reviewing some important facts about Lie groups and Lie algebras, which will be used 
in the remaining part of the paper. The exposition is very brief and follows |^6|. The reader is 
advised to refer to El 0, 0] for more details and a rigorous treatment on the subject. 



2 Cartan Decomposition and Riemannian Symmetric 

Spaces 



We will assume that the reader is familiar with the basic facts about Lie groups and homogeneous 
spaces pj. 

Notation 1 Throughout the paper, G will denote a compact semi-simple Lie group and e its identity 
element (we use I to denote the identity matrix when working with the matrix representation of the 
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group). As is well known there is a naturally defined bi- invariant metric on G, given by the Killing 
form. We denote this bi-invariant metric by <,>g- We will use K to denote a compact closed 
subgroup of G. Let L(G) be the Lie algebra of right invariant vector fields on G and similarly L(K) 
the Lie algebra of right invariant vector fields on K. There is a one to one correspondence between 
these vector fields and the tangent spaces T e (G) and T e (K), which we denote by g and 6 respectively. 
There is a direct sum decomposition g = m © t such that m = £ with respect to the metric. 

To fix ideas, let G = SU (n) and g = su(n) be its associated Lie algebra ofnxn traceless skew- 
Hermitian matrices . Then < A, B >g— trace(A' i B), A, B e su(n) (which is proportional to the 
Killing metric) represents a bi-invariant metric on SU(n). 

Definition 1 (Cartan decomposition of g) Let g be a real semi-simple Lie algebra and let the 
decomposition g = m © t, m = t 1 - satisfy the commutation relations 

[M] c e (2) 

[m,6] = m (3) 
[m,m] c 6. (4) 

We will refer to this decomposition as a Cartan decomposition of g. The pair (g, t) will be called an 
orthogonal symmetric Lie algebra pair pf], B] . 



Definition 2 (Adjoint orbit)The Lie group G acts on its Lie algebra g by conjugation Ado ■ 2 — > g 
(called the adjoint action). This is defined as follows. Given U GG, leg, then 

Mu(x) _ 4_U^pU^ u=o 

We use the notation 

Ad K (X) = |J Ad k (X). 
Adx(X) is called the adjoint orbit of X. 

For example if U S SU(ri), A £ su(n), its associated Lie algebra ofnxn traceless skew-Hermitian 
matrices . Then Adu(A) = UAU T , such that 

cxp(Adu{A)) = J7exp(A)f7 T = Ad v (exp{A)) . 

Definition 3 (Cartan subalgebra) Consider the semi-simple Lie algebra g and its Cartan de- 
composition g = m © 6. If h is a subalgebra of g contained in m, then f) is Abelian because [m, m] G t. 
A maximal Abelian subalgebra contained in m is called a Cartan subalgebra of the pair (g, 6). 

It is well known |^| that: 

Theorem 1 If f) and h' are two maximal Abelian subalgebras contained in m, then 
1. There is an element k S K such that Adk(t)) = f)'. 
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2. m = \J keK Ad k (t>). 



Remark 1 It is also well known that the homogeneous coset space G/K = {KU : U G G} admits 
the structure of a diffcrentiable manifold [[!]. Let n : G — ► G/K denote the natural projection 
map. Define o £ G/K by o — 7r(e). Given the decomposition g — m © t, the tangent space plane 
T {G/K) can be then identified with the vector subspace m. If g = m © 6 is a Cartan decomposition 
then the homogeneous space G/K = exp(m), and is called a globally Riemannian symmetric space 
||. From above stated theorem [j], the maximal Abelian subalgebras of m are all Adx conjugate 
and in particular they have the same dimension. The dimension is called the rank of the globally 
Riemannian symmetric space G/K. 



Theorem 2 Given the semi-simple Lie algebra g and its cartan decomposition g = m © t, let f) be 
a Cartan subalgebra of the pair (g, 6) and define A — exp(f)) C G. Then G = KAK. 

Proof: G = KP, where P = exp(m) = cxp(\J keK Ad k (t>)) = \J keK A<4(exp(f))) = \J keK Ad k {A) C 
KAK. Now G = KKAK = KAK. Q.E.D. 



Definition 4 (Cartan decomposition of G) The decomposition G = KAK of the semi-simple 
Lie group G = KAK , will be our most important tool in this paper. We will call this decomposition 
the Cartan decomposition of the Lie group G H 



Definition 5 (Weyl orbit) Given the Cartan decomposition g = mffit, let f) be a Cartan subalgebra 
of (g, i) containing X. We use the notation W(X) — f) p| Adx (X) to denote the maximal commuting 
set contained in the adjoint orbit of X. The set W(X) is called the Weyl orbit of X. 



3 Product Operator Basis 



The Lie Group G which we will be most interested in is SU (2 ra ), the special unitary group describing 
the evolution of n interacting spin ^ particles ( Please note that we focus on SU(2 n ) instead of U (2 n ) 
because a global phase is not of interest to us). The Lie algebra su(2") is a 4™ — 1 dimensional space 
of traceless n x n skew-Hermitian matrices. The orthonormal bases which we will use for this space 
is expressed as tensor products of Pauli spin matrices |l^] (product operator bases) . We choose to 
work in these bases because of their widespread use in the NMR literature and our desire to look at 
the implementations of NMR quantum computers. Recall the Pauli spin matrices I x , I y , I z defined 

by 
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are the generators of the rotation in the two dimensional Hilbcrt space and basis for the Lie algebra 
of traceless skew-Hermitian matrices su(2). They obey the well known relations 



where 



ll=I? = /? = -! 



1 
1 



(5) 
(6) 



Notation 2 The orthogonal basis {iB s }, for su(2 n ) take the form 

n 

Bs = 2 q - 1 Y[(h a ) ak % 

a = x,y,or z and 



k=l 



(7) 



(8) 



where I a the Pauli matrix appears in the above expression only at the k th position, and 1 the two 
dimensional identity matrix, appears everywhere except at the k th position. a^ s is 1 in q of the 
indices and in the remaining. Note that q > 1 as q = corresponds to the identity matrix and is 
not a part of the algebra. 



Example 1 As an example for n — 2 the basis for su(4) takes the form 

9=1 i{hx, Ily, hz: hx, h y , hz} 



iihxhx, hxhy, hxhz 

Ilyhx i Ilyhy , Ilyhz 

hzhx, hzhy, hzhz-} 



Remark 2 It is very important to note that the expression depends on the dimension n. To 
illustrate what this means, the expression for I 2z for n = 2 and n = 3 is 1 ® I z , and 1 ® I z ® 1 
respectively. Also observe that these operators are only normalized for n = 2 as 



tr{B r B s ) = 5 rs 2 n - 2 
To fix ideas, lets compute one of these operators explicitly for n = 2 



(9) 



which takes the form 



hz 



1 


' 1 







' 1 " 


2 





-1 




1 




" 1 
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-1 
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4 The Two Qubit Example 



Before we consider the most general case of n qubits, let us make concrete the ideas developed in 
the previous section with the help of an example. 

Example 2 Suppose we have two heteronuclear interacting qubits and the interaction between 
them produces a unitary transformation of the form 

exp(— iaJhzhz), acR. 

This is a typical scenario of two nuclear spins coupled by a scalar J coupling. Furthermore assume we 
can individually excite each spin, i.e. perform one qubit operations. The goal now is to produce any 
arbitrary unitary transformation U 6 5/7(4), from this specified coupling and one qubit operations. 
This structure appears often in the NMR situation. The unitary propogator U, describing the 
evolution of the system in a suitable rotating frame is described by 

4 

il =-i{H d + ^UiHi)U, U(0) = I (10) 

i=l 

where 

H d = 2ttJIi z I 2z 
Hi = 2irl lx 
H 2 = 2nl ly 
H 3 = 2ttI 2x 

Hi — 2TTl 2yi 

The symbol J represents the strength of the scalar coupling between I\ and I 2 . Observe that 
the subgroup K generated by {Hi}j =1 is SU{2) (g> SU{2). Therefore the unitary transformations 
belonging to SU(2) <E> SU(2) can be produced very fast by hard pulses that excite each of the spins 
individually. 

The Lie algebra 9 = su(4), has the direct sum decomposition 9 = m © t. Where 



t = span i{h Xl hyJi Zl I 2xi I 2 y,I 2 z} 
m = span i{h x I 2x , h x I 2y , h x I 2z 

hyhx, hyhy, hyhz 
Ilzhx, hzhy, hzhz} 

Please note that span in above equations and in the rest of the paper denotes all linear combinations 
with real coefficients. Using the well known commutation relations 

[A ® B, C <g> D] = [A, C] <g> (B.D) + (C.A) ® [B, D], 

and equations (|| ^), it is easily verified 

[t, 6] c £, [m, {] = m, [m, m] c t 

Therefore the decomposition 9 = m © t is a Cartan decomposition of su(4). As the subalgebra 
t = su(2) ©su(2), generates the group K — SU(2) ® 517(2), the coset space 

5/7(4) 



5/7(2)0 5/7(2) 
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is a Riemannian symmetric space. Note that the Abelian subalgebra t) generated by 

is contained in m and is maximal Abelian and hence a Cartan subalgebra of the Symmetric Space 
5[/(2)®S[7(2) ■ Therefore using the corollary [l] any Uf & SU (4) can be decomposed as 

Uf = Ki exp{-i(aihxhx + a 2 h y I 2y + a 3 h z I 2z ))K 2 . 
where K u K 2 e SU(2) SU{2). 

Now lets see how this decomposition makes obvious the choice of pulse sequences for producing this 
propogator. Note that for K y = exp(— i^h y ) exp(— i^I 2y ) , we have 

K y e^p(-iaih z I 2z )Ky 1 = exp(-iatili x l 2x ) 

and similarly for K x — exp(— i%Ii x ) exp(— i%I 2x ) we have 

K' 1 exp(-ia 2 Ii z I 2z )K x = c^p(~ia 2 h y I 2y ). 

This makes transparent, as to how one should generate the unitary transformation Up above using 
the unitary evolution in equation (|l0|). This can be summarized by writing Up as 

Up = Kx ■ K v ■ exp(-iaili z l 2z ) ■ Ky 1 ■ K^ 1 ■ exp(-ia 2 h z I 2z ) ■ K x ■ exp(-ia 3 Ii z I 2z )) ■ K 2 . (11) 

The unitary propogators K Xl K yi K\ and K 2 can be produced by selective hard pulses, i.e. one qubit 
gates. 



Remark 3 In p7| , it was shown that synthesizing Uf using the decomposition given above, is 
indeed the fastest way to generate Up. Time optimality is an important consideration in presence 
of decoherence. Observe that the Hamiltonians l\ x l 2x and I\ y I 2y are the elements of Weyl orbit of 
Iizhz under the adjoint action of SU(2) <g> SU{2). 



Now lets proceed to the general case of n qubits. 



5 How does [SU(2)] m sit inside SU{2 n ) 



In this section we would like to use the machinery developed in the previous section to decompose 
the group SU(2 n ), into smaller unitary transformations, which can be produced. This will be 
achieved through successive Cartan decompositions of SU(2 n ) into smaller and smaller unitary 
transformations till we are only left with unitary operations corresponding to one and two bit 
operations. 

Main Idea: Suppose we are given n qubits. In NMR quantum computing this will be a network of n 
coupled spin \ particles. Let us label these qubits as qi, q 2 . . . q n . Let Hi denote the two dimensional 
Hilbert space of a single qubit. Similarly let H n denote the 2™ dimensional Hilbert space of n qubits. 
These are related by 

H„ = Hi ® Hi ® . . .Hi. 
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To understand the decomposition of SU(2 n ), suppose we can generate an arbitrary unitary trans- 
formation on the Hilbcrt space H n _i of qubits q\, q<i . . . q n -i i-c an arbitrary element of SU(2 n ~ 1 ). 
Also assume, we can independently manipulate the n th qubit. Furthermore assume we have a two 
qubit gate that will act on the qubit q n and q n -\. In the context of NMR, this corresponds to evo- 
lution under a coupling between the spin q n and q n -i- To be more specific let this unitary evolution 
be of the form 

exp(-iaJI {n _ 1)z I nz ), 

caused due to a scalar J coupling. 

We now explicitly build any U £ SU(2 n ) from elements of SU(2 n ' 1 ) ®)SU(2) and this interaction 
between q n -i and q n . Notice this will essentially solve our main problem because to generate 
any element V £ SU(2 n ^ 1 ) we will use this divide and conquer strategy again, building V from 
SU(2 n ~ 2 ) ® SU(2) and interaction between (n — l) th and (n — 2) th qubit and so on. Let us look at 
the geometry of the situation 

Notation 3 The Lie algebra su(2") consists of the elements 

su(2") = span {A® I x , B ® I y ,C <g> I z , D ® 1 , il nx , il ny , U nz \A, B, C, D e su(2™ -1 )} 
The Lie algebra su(2 n ) has a direct sum decomposition into the following two vector spaces 



We emphasize again that span denotes the vector space generated by the elements of the set by 
taking linear combinations over the field of reals. 



Lemma 1 The vector space SU{(2") is a Lie algebra of dimension 2x4™ 1 — 1 such that 



Proof: The proof is a direct consequence of following commutation relations. Let A, B, C £ 



su { (2") = span {A® I Z ,B® 1 ,U nz \A,B e su(2"" 1 )}, 
su m (2") = span {A® I X ,B <g> I y ,i \A,B esu(2™- 1 )}. 



exp(su { (2™)) = SU{2 n - r ) ® SU{2 n ' 1 ) ® U(l) 



su(2 



if A^B 



[A®I Z ,B®I Z ] = 
[A®I Z ,B®1] = 
[A ® 1 , B <8> 1 ] = 



C<8>1 

C®I Z 

C®1 



[A®I z ,I nz ] = 



0; 
0; 
0. 



[A ® 1 , I nz ] = 
[A®I Z ,A®1\ = 



Now observe that I nz and elements of the type 




1 










and 











1 



commute. Therefore the Lie Algebra 



su { (2") = su(2™- 1 ) © su(2™- 1 ) © u(l) 



and the last part of the lemma follows. 



Q.E.D. 
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Theorem 3 The decomposition su(2 n ) = su m (2 n ) © SU{(2") is a Cartan decomposition of the Lie 
algebra su(2 n ) 

Proof: From equation (||), it is clear that SU{(2™) _L su m (2"). We have already shown in Lemma |l| 
that su{(2 n ) is a Lie subalgebra of Su(2 n ). All that we need to show is that the commutation rules 

[ S u m (2"),su { (2")] = 5u m (2") 

[5U m (2"),5U m (2")] C 5U 6 (2"). 

are satisfied. The proof follows from the following commutation relations and the fact that su(n) 
is semi-simple (implies [su(n) , su(n)] — su(n)). Let A,B,C G su(2™ _1 ), c G R and a G {x,y}. If 
A ^ B, then 



If A = B then 



[A®I a ,B$ 




= C7®1 ; 


[A®I x ,B<k 


)Iy] 


= C®I Z ; 


[A®I X ,B§ 


5 4] 


= C®I y ; 


[A®Iy,B$ 


5 4] 


= C®I X ; 


[A®I a ,B<k 


51 ] 


= C®I a . 


[A® I x , Ac 


^Iy] 


c il nzi 


[A®I x ,At 


S>4] 


C Uny, 


[A®Iy, AC 


3>4] 


c ilnxi 


[A® I a ,Ai 


2)1 ] 


= 0. 



Corollary 1 Any C/ G SU(2 n ) has a decomposition 

where K 1 ,K 2 G exp(su e (2")) = 5t/(2"- 1 ) ® SUIT 1 - 1 ) <g> C/(l) and A = exp(Y), for some Y € J), a 
maximal Cartan subalgebra of the pair (su(2™), sut(2™)). 

Proof: We have already shown that su(2") = siif(2™) ©su m (2") is a Cartan decomposition. Hence 
the proof follows from theorem |^. Q.E.D. 

The key observation is that K\, K2 G exp(su{(2 ra )) again has a decomposition. This is stated in the 
following theorem. 



Notation 4 Let us denote 

su n (2 n ) = span {A® I Z \A G su(2™ -1 )} 
su {0 (2") = span {A <8> 1 \A G su(2"- 1 )} 
su { (2") = span {A <g 1 , A® J»| A G su(2"~ 1 )} 

Notice SU{(2") is just SU{(2") © u(l) (I„ 2 commutes with eveything in SU{(2") and generates u(l)) 
and therefore 

exp(su t (2™)) = exp(su { (2")) ® Z7(l). 
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Theorem 4 The direct sum decomposition sut(2") = sUfco(2™)©sUfci(2™) is a Cartan decomposition 
of the Lie algebra SU{(2"). 

Proof: Recall that 

su { (2") = span {A ® I z , B ® 1 |A, fl 6 su(2™" 1 )} 
If 4,B,Ce su(2"- 1 ) such that A ^ S, then 

[>1(8»J 2 ,B(8»1] = C®I Z ; 
[A®I Z ,A®1] = 0; 
[A<g> J z ,B(g)I z ] = C®1 

Therefore from the above relations and the fact that su(n) is semi-simple, it can be verified that 

[suei(2"),BUn(2")] c su t0 (2") 
[su ei (2"), S u {0 (2")] = su (1 (2") 
[su eo (2"), S u {0 (2")] c su (0 (2") 

Hence the result follows. Q.E.D 

Corollary 2 Any U E cxp(suf(2 n )) has a unique decomposition 

U = K X AK 2 , 

where K\,K 2 E SU(2 n ^ 1 ) and A — exp(Y), for some Y £ f, a Cartan subalgebra of the pair 
Proof: The proof follows directly from theorem ||. 

The above stated theorems, therefore give us a recursive decomposition procedure for an element in 
SU(2 n ). To summarize what we have accomplished till now. 

Corollary 3 Any U E SU(2 n ) has a decomposition 

U = Ki- exp(Zi) • K 2 • cxp(r) • K 3 • exp(Z 2 ) • K 4 , 

where K x , K 2 , #3 , #4 € St/^™" 1 ) <g) U(l), Y € f), a Cartan subalgebra of the pair (su(2"), su t (2™)) 
and Zi, Z 2 G f, a Cartan subalgebra of the pair ^sU{(2"),sueo(2")^ . 

Proof: The result follows directly from corollaries [l] and |2|. 

Remark 4 We know how to produce unitary transformations K\, K 2l K3, K4 as these belong to the 
group SU{2 n - 1 ) O 17(1), a subgroup of SU(2 n - 1 ) ® SU{2). All that needs to be shown now is that 
we can produce unitary transformations of the form exp(Zi), exp(^2) and exp(Y), where Y £ t), 
a Cartan subalgebra of the pair (su(2 ra ), sU{(2 ra )) and Z%,Z 2 E f, a Cartan subalgebra of the pair 

(sut (2"), sUfo(2")^ .This will be achieved by using the coupling between the n th spin and the network 

of coupled Ti—l spins. We will show that the adjoint action of SU(2 n ~ 1 ) eg) SU(2) on the coupling 
Hamiltonian (I( n —x\ z I n z) contains commuting elements whose span is the whole Cartan subalgebra 
f). Same applies for the Cartan subalgebra f (Recall we used this fact in the two spin case, where we 
generated any element of the Cartan subalgebra, of the form axli x l 2x + a. 2 l\ y l 2: u + ctshzhz by the 
adjoint action of SU{2) ® SU{2) on I\ z I 2z ). Now using the result of corollary (0) we will be able to 
generate any unitary transformation on n qubits. 
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5.1 Generating Cartan subalgebras from Weyl orbits 



Let us now characterize a Cartan subalgebra for the pair (su(2"),su E (2")) and (su t (2™), su {0 (2")J , 

in the our product operator basis. This will be done in a recursive way. For this we use the following 
notation 

Notation 5 Let 

a(2) = ±{ih x l2x,ihyh y ,ihzl2z} 
«(2) = o(2) 

a(n) = {±U nx ,A®I x \A e s(n- 1)} 
b(n) = {A® I Z \A e s(n - 1)} 
s(n) = {s(n - 1) (g> 1 , a(n)} 
F)(n) = span a(n) 
f(n) = span b(n) 



Remark 5 To avoid confusion, we remind the reader that for A e s(n — 1), the elements A ® 
and A ■ I nx represent the same object. 

Example 3 To fix ideas, we give here the explicit expressions for a(n) and b(n) for n — 3 and n = 4 
a(3) = ± i{h x , Ii x hxh x , Iiyhyhx, hzhzhx} 

tt(4) = ± i{Ii Xl h x lAxTllxhxIi X ThyhylAxThzhzlAxTh X hxh X li X Thyhyl3xIi X Tllzl2zIs X Ii X } 
b(3) = ± l{Il x l2 X l3z, hyl2yl3z, hzjlzlzz} 

b(4) = ± i{h x Iiz, Ilxhxliz, Ilyhyliz, Ilzhzliz, Ilxl2 X h X Iiz, IlyhyhxI^z, IlzhzhxI^z} 

Lemma 2 The subspaces \) (n) and f (n) are a Cartan subalgebra for the Lie algebra pair (su(2™ ) , SU{ (2 
and ^su{(2™),su{o(2")^ respectively. 

Proof: We first show that t)(n) is a Cartan subalgebra of (su(2™), SU{(2™)). The proof is inductive. 
We saw that f}(2) is a Cartan subalgebra of the pair (su(4), su(2) © su(2)). We assume that the span 
of s(n— 1) is maximally Abelian in su(2™ _1 ) (easily verified for s(2)). We now show l)(n) is Abelian. 
Observe, it suffices to prove that all the elements of a(n) commute. The commutation relations 

[I nx ,A-I nx ] =0, Aes(n-l), 

[B ■ I nx , A ■ I nx ] = 0, A,BGs(n-l), 

and 

[B, A ■ I nx ] = 0, A,Bes(n-l), 
suffice to show that a(n) and s(n) are Abelian too. 

To show t)(n) is maximally Abelian in su m (2 n ), we use induction again. Recall that 
su m (2 n ) = span{A ® I X ,A® I y ,iI nx ,iI ny \A esu(2™" 1 )}, 
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and I nx does not commute with I ny , I nz and A® I y , A® I z , where A 6 su(2" _1 ). The set 
s(n — l).I nx is maximally Abelian in su(2™ _1 ) <E> I x (Because s(n — 1) is maximally Abclian in 
su(2 n_1 ) by assumption). Therefore the span of a(n) = {±iI nx ,A ■ I nx \A e s(n — 1)} is maximally 
Abclian in su m (2") and the span of s(n) = {a(n),s(n — 1) (g 1 } is maximally Abelian in su(2 n ). 
Hence the induction argument is complete. 

The proof f(n) is a Cartan subalgebra of the pair (sue(2™), sU{o(2 n )^ follows exactly on same lines, 
and we leave it for the reader to verify it. Q.E.D. 

Remark 6 We now show that we can generate any element of the Cartan Subalgebras f)(n) and 
f(n), via the adjoint action of SU(2 n ~ 1 ) (g SU(2) on the coupling Hamiltonian I( n -i)zlnz- To be 
more precise, any element Y G f) can be written as 

p 

Y = Y,frAd Kj {I (n _ 1)z I nz ), (i, > 0, Kj e SU(2 n - 1 ) ® 5(7(2), 
i=i 

such that AdKjil^n-^zInz) all commute (elements of a Weyl orbit under the adjoint action of 
SUIT 1 - 1 ) ® SU(2) ) and therefore 

exp(y) = iF^Kj ■ eMhhn-DzInz) ■ K). 

Lemma 3 Let K = SU(2 n ^ 1 ) (g) SU(2), then the adjoint orbit AdK{I( n -i)zInz) contains the sets 
s(n — 1) (g) I x and s(n — 1) ® I z = b(n). 

Proof: The proof is again inductive. Note by definition 

s(n - 1) <g) I z = {s(n - 2) ® 1 ® 7 Z , s(n - 2) O 7 X ® 7 Z , ±i/ (n _ 1)a ./ nz }. 

We assume that the statement of the lemma is true for n — 1, i.e., if = S'?7(2" _2 ) <g> 5C/(2), then 
Ac?h (/( n -2)z7(n-i)z) contains the set s(n— 2)(g7 x . Therefore the adjoint orbit AdK{I(n-2)zI(n-i)zInz) 
contains the set s(n — 2) <g> 7 X <g) 7^ (Apply i7 to I( n -2)zl(n-i)z an d don't do anything to last spin 
(qubit)!). Now observe for 

Ui = exp(-i7r7(„_ 2 ) z 7( ri _i) !/ ) 

and 

U 2 = exp(i-7 (n _ 1) j / ) 

both belonging to if, we have 

U2U\{I(n-\)zInz)U\u\ = 7( n _ 2 ) z 7( n _i) z 7 nz . 

Since the term ±7(„_ 1 ) x 7 nz are present in adjoint orbit ^4<iA:(7(„_i) z 7„ z ), (By selective | rotation 
of n — 1 qubit) we deduce that the adjoint orbit Adx(7(„_i) z 7„ z ) contains the set 

o(n - 1) ® I z = {s(n - 2) ® 7 X (g) 7 Z , ±7 (n _ 1)x 7 nz }. 

Finally observe terms of the form s(n — 2) (g 1 ®/ z can be produced by adjoint action of K\ = 
exp(— iirl{ n -\) x lny) on the set s(n — 2)(g7 x (g7 x , which produces Hamiltonians of the form s(n — 2) ® 
1 ® 7 Z . Thus we can generate the whole set s(n — 1) ® 7 Z . That we can also generate s(n — 1) ® 7 X , 
is also obvious (by selective \ rotation of the n th spin). Q.E.D. 
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Remark 7 Note that the Hamiltonian ±I nx can be produced by selective excitation of the n th 
qubit. Thus we can produce any Hamiltonian from the set {s(n — 1) ® I x , ±I nx } = a(n). We have 
therefore shown that we can produce any Hamiltonian Y and Z belonging to the set a(n) and b(n) 
and as there positive span is the Cartan subalgebra f)(n) and b(n) we can generate any Hamiltonian 
belonging to \)(n) and b(n) by the adjoint action of the group SU(2 n ~ 2 ) SU(2) on the coupling 
Hamiltonian and we are done. We complete the inductive proof with the following lemma. 



Lemma 4 Let K = SU(4)®SU(2), then the adjoint orbit Adxihzhz) contains the sets 0(2)0/3; 
s(2) <8> I x and a(2) I z = b(3). 



Proof: Recall 



d(3) — ± i{I 3x , Ilxhxhx, hyllyl^x, Ilzhzhx} 

b(3) = ± i{hxhxhz, hyhyhz, hzhzhz} 



The result then follows from the following relations 

ey^{-iTTl lx I 2 y)e^{-i'Kal2zhz)^>{i^hxhy) = cxp(-iTTaIi x I 2 xhz) 
exp(-i^iu) cxp(-i^I 2z ) expi-iirahxhxhx) exp(i^I 2z ) exp(i^7 u ) = cxp(-z7ra/i !y /2j / ^) 

7T 7T 7T 7T 

cxp(-2 - /ij,) exp(-i-I 2y ) exp(-ianli x l 2x l 3x ) exp(i-I 2y ) exp(i-I ly ) = cxp(-iTraI lz I 2z I 3z ) 

which shows how to generate b(3). Now terms in o(2) I x can be generated by selectively rotating 
the third spin. Q.E.D. 

To review the ideas developed and to conclude the paper, let us explicitly work out the synthesis 
of any unitary propogator in a linearly scalar coupled three spin system. 



Example 4 Consider three heteronuclear spin | nuclei, coupled through scalar coupling . One has 
the ability to selectively excite each of the nuclei. Thus the system evolution in a suitable rotating 
frame is approximated by 

6 

U=-i{H d + Y,u l H l )U, (12) 
i=i 

where 



H d 


= 2irJ {h z hz + hzhz) 


Hi 


= 2ttIix 


H 2 


- 27!-/^ 


H 3 


= 2ttI 2x 


H A 


= 2irl 2y 


H 5 


= 2irl 3 x 


H 6 


= 2-Kl 3y 



Observe that the subgroup generated by is SU(2) <8> SU(2) SU(2). Thus any unitary 

propogator U can be decomposed as 
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U = Ki exp(-i(f3ih x I 2 xhx + (hhyhyhx + Pzhzhzhx + Pihx))K 2 

where 

Ki = Pi exp(-i(ailixhxhz + ot2h y h y hz + azhzhzhz))P2 

and 

K 2 = P3exjp{-i("fih x hxhz + ^hyhyhz + ~i3hzhzhz))Pi, 

where Pi , P 2 , P3 , P4 belong to SU(4) <g> SU(2), with 5[/(4) denoting arbitrary transformation on 
spin 1 and 2 and SU(2) representing selective local transformations on the spin I3. 

Any element of SU (4) can be produced by decoupling the spin I3 (rapidly flipping the spin j^] ) 
and then using the coupling I\ z I 2z as in example |[ The Hamiltonians belonging to the Cartan 
subalgebra f)(3) and f(3) can then be produced by using the coupling I 2z Hz as illustrated in lemma 



Remark 8 We have therefore demonstrated constructive controllability in a network of coupled 
spins. The crucial thing to note is that we have dealt with the worst case scenario. The case we 
have treated here is of a network of linearly coupled spins. There was no direct coupling between 
the n th spin and say k th spin for k < n (In terms of quantum computing, if k < n — 1, we are 
not allowed to let k th and n th qubit interact with a two qubit gate directly). However interactions 
can be mediated through the other spins and this is what our constructive procedure is doing. It is 
also not difficult to see now that as long as we have a network of spins which is connected we have 
constructive controllability. 



6 Conclusions and Future Work 



Our main goal in this paper has been to put the design of quantum computers and pulse sequences 
in NMR on a geometrical footing. We have produced here a parameterization of the group SU(2 n ) 
in terms of unitary transformations produced by one and two qubit gates using successive Cartan 
decompositions of SU(2 n ). The first issue we would like to draw attention to is time optimality. 
We would like to emphasize that our constructive procedure of producing an arbitrary unitary 
transformation is only time optimal for the two spin case . Recall that in the course of synthesizing 
a propogator almost all of the time is spent during the interaction of qubits and this corresponds 
to the two qubit gates or evolution of couplings during NMR. The one qubit operations can be 
produced relatively fast by external selective hard pulses. Thus from a practical viewpoint, it is of 
utmost importance that not only do we have a constructive procedure for synthesizing an arbitrary 
unitary transformation, but one which is time optimal to minimize the effects of decoherence which 
are always present. In p7| , we developed the theory for time optimal control of spin systems and 
computed the time optimal control laws for two spin system (as in example^) for any kind of coupling 
between the spins. The two spin case is special and elegant because su^]®SU(2) na PP ens to be a 
Riemannian symmetric space. The results of time optimality for two spin systems do not extend in 
a natural way to higher spin systems as the coset space ^j^^)]®" does n °t have a symmetric space 
structure. Finding time optimal control laws for spin networks with more that two spins will require 
and inspire further developments in geometric control theory. 

Another question of imminent interest is - how to transform between different paramctcrizations of 
SU(2 n ). Suppose U — SU{2) , then we can express any element of SU(2) in the following two ways 

U = exp(-iail x + a 2 I y + a 3 I z ), ai , a 2 , a 3 <G R (13) 
U = exp(-i!3 1 Ix)e X p(-iP2ly)exp(-0 3 Ix), ft,ft,ftel (14) 
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It is well known, how to transform between these and other ways of expressing an element of SU (2). 
It will be very interesting to find the transformation that takes our description of SU(2 n ), which is 
in the same spirits as the last of the above descriptions, to the other standard parameterizations of 
SU(2 n ). 
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